Coupled dynamics of electrons and phonons in metallic nanotubes: 
current saturation from hot phonons generation 
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We show that the self-consistent dynamics of both phonons and electrons is the necessary ingre- 
dient for the reliable description of the hot phonons generation during electron transport in metallic 
single- wall carbon nanotubes (SWNTs). We solve the coupled Boltzmann transport equations to 
determine in a consistent way the current vs. voltage (IV) curve and the phonon occupation in 
metallic SWNTs which are lying on a substrate. We find a good agreement with measured IV 
curves and we determine an optical phonon occupation which corresponds to an effective tempera- 
ture of several thousands K (hot phonons) , for the voltages typically used in experiments. We show 
that the high-bias resistivity strongly depends on the optical phonon thermalization time. This im- 
plies that a drastic improvement of metallic nanotubes performances can be achieved by increasing 
the coupling of the optical phonons with a thermalization source. 

PACS numbers: 61.46.Fg, 63.20.Kr, 65.80.+n, 72.15.Lh, 73.63.Fg 
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I. INTRODUCTION 

Much interest is currently devoted to the study of elec- 
tronic transport properties in carbon nanotubes (CNTs), 
both metallic and semiconducting. The main reason 
is the possibility to use CNTs in electronic integrated 
circuits thanks to catalizator-assisted on-site growth^. 
Semiconducting CNTs are envisageables as new compo- 
nents for transistors and logic circuits 3 ' 4 . On the other 
hand, metallic CNTs are particularly suited as a new 
type of interconnects due to their small dimensions and 
the large electron current density (^10 9 A/cm 2 ) they can 
support. The current vs. voltage (IV) curve of metallic 
single-wall carbon nanotube (SWNTs) has been recently 
measured by several groups^J. For voltages > 0.2 V, 
they observe a sudden increase of the resistivity which 
is due to the scattering with optical phonons. In long 
tubes, this leads to a saturation current of ~ 25 //A for 
voltages > 5 V. Such behavior limits the performances 
of metallic SWNTs as interconnects. The understand- 
ing of this phenomenon is, thus, a crucial step toward 
finding methods to boost SWNTs performances and has 
important technological consequences. 

Two recent papers^ suggested the possibility that, at 
high bias, the electron transport induces an anomalously 
high optical-phonons occupation (hot phonons) which, 
in turn, induces an increase of the resistivity. In Ref. |g, 
this hypothesis is formulated on the basis of the com- 
parison between scattering lengths (obtained from IV 
measurements Ml) and electron-phonon coupling val- 
ues (obtained from ab-initio calculations^ and inelastic 
X-ray scattering measurements^). However, the extent 
of Ref. conclusions is limited, because Ref. does not 
provide neither a scheme to reproduce experimental IV 
curves, nor a reliable quantitative determination of the 
phonon occupation. The conclusions of Ref. apply to 
tubes which are lying on a substrate. This is the very 
typical situation encountered in experiments. 

Ref. 9 reports the comparison between IV curves mea- 



sured on SWNTs in two different situations: when the 
tubes are lying on a substrate and when the tubes are 
suspended between the two electrodes. Ref. |9j concludes 
that hot phonons are present in the tubes which are sus- 
pended but are absent in those which are lying on a sub- 
strate (for bias < 1 V) . The presence of hot phonons is 
clearly demonstrated by the experimental observation of 
a negative differential resistance in suspended tubes. On 
the other hand, the absence of hot-phonons is inferred 
from a simplified theoretical model in which the electron 
scattering length (at zero temperature) is a fitting pa- 
rameter. The claim of absence of hot phonons, for the 
tubes lying on a substrate, is contradicting the conclu- 
sions of Ref. |8|, where the scattering length is obtained 
from ab-initio calculations of the electron-phonon cou- 
pling. 

Concluding, a quantitative model to describe the 
anomalous phonon heating is not available. More over, 
given the contradicting conclusions of Refs. Islil it is 
not clear whether hot phonons generation is present for 
tubes lying on a substrate. This last point is particularly 
relevant. In fact, if the high-bias resistance measured for 
tubes lying on a substrat o 5 i 6 i 7 is due to hot phonons — , 
one could conceive some methods to reduce the optical 
phonon temperature and thus diminish the high bias re- 
sistance. This would not be possible if hot phonons were 
not present in tubes lying on a substrate 9 -. 

In this paper, we solve the coupled Boltzmann trans- 
port equations for both phonons and electrons to deter- 
mine in a consistent way the IV curve and the phonon 
occupation in metallic SWNTs. The quantitative deter- 
mination of the anomalous phonon heating during elec- 
tron transport allows us to settle the debate on the pres- 
ence/absence of hot phonons in nanotubes lying on a 
substrate. We are interested in the high bias (> 0.2 V) 
region, where the transport properties are determined by 
the scattering with optical phonons. We use the zero- 
temperature electron scattering lengths obtained from 
ab-initio calculations based on the density functional the- 
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FIG. 1: Scheme of the model used in this paper (see the text). 



times. T^p and are phonon scattering times, i.e. the 
average time a phonon lives before emitting an electron- 
hole pair. For metallic SWNTs, it can be shownSS that 
T ep = T 6s/4 an d t% = T ^/2- The optical phonons can, 
also, scatter into acoustic phonons with a characteris- 
tic thermalization time Tth ~ 5 ps. Both the acoustic 
phonons and the substrate are thermalized at room tem- 
perature, and are acting as a thermal bath. The present 
model is not valid if the tube is suspended between the 
two electrodes, as in some of the experiments of Ref. U 
In this case, the acoustic phonons are not thermalized 
with the environment and their occupation should be de- 
termined self-consistently. We are not considering this 



ory 



8.10.12 



II. DESCRIPTION OF THE MODEL 
A. GENERAL CONSIDERATIONS 

In a metallic SWNT, there are only two phonons which 
can be generated by transport electrons via a back- 
scattering process 8 . Using the graphene notation, they 
correspond to the F,2 g LO (here LO stands for longitudi- 
nal along the tube axis) mode at T and to the mode at 
K (equivalent of the A^ at K')i^. For simplicity, we label 
them r and K. r£ and t^ are the corresponding elec- 
tron scattering times, i.e. the average time an electron 
can travel before emitting a phonon. t& s = lb s /vp, where 
lbs is the scattering length and vf = 8.39 x 10 7 cm/s 
is the Fermi velocity^. According to precise ab-initio 
calculations^, at zero temperature r£ = 538 and t£ = 
219 fs^ (for tubes with a diameter of 2.0 nm, which is 
the diameter typically used in experiments &£il) . 

Electron transport can induce not-thermal high- 
occupation of phonons if the generation time of the 
phonons (t£ s and r^) is shorter than their thermaliza- 
tion time. For metallic SWNTs, the optical-phonons 
thermalization time Tth is due to anharmonic scatter- 
ing into acoustic phonons — . Time-resolved terahertz 
spectroscopy measurements^, on graphite, reported an 
optical-phonon thermalization time of 7 ps. This value 
is comparable with the optical-phonon lincwidth of di- 
amond (the other carbon polymorph) ~1 cm" 1 , which 
corresponds to a relaxation time of 5 ps, and is due to 
phonon-phonon scattering^ 7 -. A Tth value of the same 
order is also obtained by preliminary calculations of 
phonon-phonon scattering on graphite and na notubes — 
done with the ab-initio methods of Refs. ll7ll9i Conclud- 
ing, we expect Tth to be of the order of 5 ps. Since tR 
and t^" are much smaller than Tth ~5 ps, hot-phonon 
generation is a priori expected to occur during electron 
transport in metallic SWNTs. 

In the following, we consider the model schematized 
in Fig. n Transport electrons can scatter with the opti- 
cal phonons T, K and K' (K' is equivalent of K and 
is not shown), and are the electron scattering 



B. IMPLEMENTATION 

To study the evolution of the optical-phonons occu- 
pation during electron transport, we consider the Boltz- 
mann transport equation. In a metallic SWNT the elec- 
tron gap is zero for two equivalent k-points corresponding 
to K and K' in the graphite Brillouin zone. Two electron 
bands are crossing at the Fermi energy at both K and 
K'. For the energy range of interest ( ± 0.5 eV around 
6f), the two bands can be considered linear with slope 
±h\vp\, where vf is the Fermi velocity. fL/R.(k,x,t) is 
the probability to have an electron coming from the left 
or right contact (L,i?) with momentum k (fc = at the 
band crossing) , at the position x. States near K or K' 
have the same distribution /. / evolves in time (i) ac- 
cording to 



[d t ± v F d x - eS/nd k }f L /R = [dtfh 



(1) 



where vp is the Fermi velocity and e (e > 0) is the elec- 
tron charge. For an applied voltage V, we consider a 
uniform electric field £ = —V/L, where L is the tube 
length. The collision term [dtf] c describes the electron 
scattering. The probability distribution is given by the 
stationary solution of Eq.^ As explained in Ref . 0, [dtf} c 
can be considered as the sum of three terms: (1) the elas- 
tic scattering [9 t /i]e = {v F /l e )UR - /l), where l e is the 
elastic scattering mean free path; (2) the forward scatter- 
ing from optical phonons [dt fiefs', (3) the back scattering 
from optical phonons [dtfh\bs- Following Refs.|7J|9j we as- 
sume l e = 1600 nm. We neglect [dtfh\js since it does 
not change the propagation direction of electrons and, 
thus, should have a minor effect on the current. Now, we 
consider an optical phonon v = I\ K, with occupation 
probability n u (k,x,t) (for the K phonon, k stands for 
K+fc) . n evolves according to 



[dt + v v ph {k)d x ]n" = [d t n" 



(2) 



where v v vh is the phonon velocity and [dtn v ] c is the col- 
sign^) 2.9 x 10 5 cm/s and 
sign(fc) 7.2 x 10 5 cm/s, from ab-initio calcula- 
tions^, where sign(fc) is the sign of k. We remark 
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that, given the small phonon velocity (with respect to 
Vp) the optical phonons do not have the time to dif- 
fuse along tubes longer than 100 nm (the thermaliza- 
tion scattering lengths are lf h = vE h T t h ~ 15 nm and 



tfk = v aL T th ~ 36 nm). This was already pointed out 
in Ref. [3 Thus, the resulting IV curve is only slightly 
affected by the exact value of \v P h\- For given x and t, 
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where fc ± = k ± (u) v /vp), fh stands for //,(&), the vari- 
ables x and t are omitted, and v = T, K. tjj v is the 
phonon pulsation, ftio T — 196.0 meV, fiu; K = 161.2 meV 
and the u v dependence on k can be neglected. Notice 
that Eq. 0] can be obtained from Eqs. 1 1121 and im- 
posing the conservation of energy and momentum in the 
back-scattering processes. We recall that, to be consis- 
tent with the Boltzmann treatment, the scattering times 
t£ s and t" of Eqs. and ^ (Fig. are computed with 
zero phonon occupation^. 

We impose the equilibrium distributions as boundary 
conditions: 



f L (k,0) = f R {-k,L) = {exp[nvFk/{k B T)] 
n u (k > 0,0) = n u (k < 0,L) = 0, 



I}" 



where T = 300 K and ks is the Boltzmann constant. 
The current is given by: 



I = 



4e f 



The stationary solution of Eqs. ^ an d El is found by 
numerical integration in timeSi. 

The present approach is similar to that of Refs. Jala 
but has two important differences. First, in Refs. |a|6| 
the phonon occupation is considered to be thermalized 
at room temperature (n ~ 0). Thus, in Refs. 00 Eq.^ is 
not taken into account and, in Eq. |3J the terms depend- 
ing on n are not considered. Second, in Refs. |5||g, the 
scattering time Th s is considered as a parameter and its 
value is fitted to reproduce the experimental IV curves, 
supposing that hot phonons are absent. In the present 
work, we do not make any assumption on the phonon 
occupation n and we determine n by solving the trans- 
port equation. Moreover, the zero-temperature scatter- 
ing times t£ s are not fitted to recover a better agreement 
with measurements, but are fixed to the values obtained 
from ab-initio calculations &i£Li£. 
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FIG. 2: (Color online) Lower panel: current vs. voltage 
characteristic of a 300 nm long carbon nanotube (CNT). 
Calculations are done fixing the phonon occupation to zero 
(cold phonons) or allowing phonons to heat up (hot phonons, 
Tth = 5.31 ps). Measurements are from RefslaU Upper panel: 
phonon effective temperature of the Raman G peak (see the 
text) as a function of the voltage. 



III. RESULTS 

A. CURRENT VS. VOLTAGE 

In Fig. El we show the comparison between the mea- 
sured IV curves and our calculations done imposing the 
phonon occupation n — (cold phonons) or allowing 
the n to vary according to Eqs. |2 El an d 0] using 
T t h = 5.31 ps (hot phonons). This value corresponds 
to an anharmonic contribution to the phonon linewidth 
(full width) 7 =1 cm -1 , which is expected from inde- 
pendent calculations and measurements ^LiLiS (see dis- 
cussion above). Imposing n = the resulting IV curve 
dramatically underestimates the experimental resistivity. 
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FIG. 3: (Color online) Current vs. voltage curve calculated 
for various tube lengths (L). Tth = 5.31 ps (see the text). 
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FIG. 4: (Color online) High-bias resisitivity for various tubes 
(see the text). Measurements are from Ref|7j. Inset: the 
scattering length l sc (defined in Eq. as a function of Tth- 
The dashed line correspond to the experimental value l 3C — 
8.5 nm (see the text). 



On the contrary, if the phonon occupation n is deter- 
mined self- consistently the IV curve nicely reproduces the 
experimental data. Already at 0.2 V the hot-phonon and 
the cold-phonon curves significatively differ. This implies 
that the hot phonon generation dominates the resistivity 
behavior already at V>0.2 V, also for tubes lying on a 
substrate. 

In Fig. 13 we show our calculations for tubes with dif- 
ferent lengths L. The dependence of the IV-cur ve on L is 
similar to that observed experimentally in Refs. I(l7l De- 
spite the presence of the hot phonons, for bias > 0.5 V 
the differential resistance dV/dl is almost independent 
from the voltage. 

We computed the dV/dl at high bias (in between 
1.0 and 1.4 V) for tubes with different lengths L 
(Fig. Beside Tth =5.31 (7 = 1 cm -1 ), we also con- 
sidered T t h =2.65, 0.885, and ps (corresponding to 
7 =2, 6, +00 cm -1 ), to investigate whether a significant 
reduction of the resistance can be achieved by a reduction 



of T t h- dV/dl is clearly linear with respect to L and is, 
thus, possible to define a linear resistivity p. For a one- 
dimentional channel with 4 sub-bands, in the incoherent 
limits : 



1 dV 
P= L-dI 



^sc 4c^ 



(5) 



where h/(2e 2 ) = 12.9 kSl is the quantum of resistance and 
l sc is the effective scattering length. Eq. [S] can be used 
to obtain l sc from the measured resistivity p. Using this 
approach, several independent experimental works 5 6 7 9 
found a scattering length l sc of the order of 10 nm, at high 
bias, for tubes lying on a substrate. We remark that l sc 
is not to be confused with the lbs of the present paper. 
While l^s is calculated for zero phonon occupation, l sc is 
an effective scattering length that includes the effects of 
the finite phonon occupation. To compare our calcula- 
tions with the measured 1 SC1 we make a linear fit of the 
data in Fig.0J we obtain p and, thus, l sc through Eq. 
This is done for the different Tth values. For Tth = 5.31 ps 
(which is our best Tth estimate) we obtain l sc =12 nm. 
Again, this compare very well with l sc =8.5 nm which 
we extracted from the experimental values^ reported in 
Fig.H 

Given the fact that the IV curves strongly depend on 
the optical phonons occupation (Fig. [21, we expect that 
the thermalization time r t h might play a crucial role. In 
fact, a smaller T t h induces a lower optical phonons occu- 
pation (Fig. . In turn, this results in lower resistivity 
or a higher scattering length l sc . Indeed, for T t h =2.65, 
0.885, and ps, we obtain l sc — 22, 42, and 108 nm, 
respectively (inset of Fig. QJ. This strong dependence 
of l sc on T t h implies that even a small decrease of T t h 
induces an important amelioration of the nanotube per- 
formances. A higher optical phonon thermalization can 
be obtained, e.g., by using a substrate characterized by 
vibrational frequencies similar to those of a CNT— . 

Finally, the calculations done with T t h — are equiv- 
alent of imposing n v — (cold phonons) in Eqs. |3 01 
For n = 0, it is usually assumed that l sc is a mea- 
sure of the optical phonon back-scattering length, i.e. 
(hc)^ 1 — Uba) -1 + Cbs) -1 - Indeed, from our calcula- 
tions for n = (r t h = 0) l sc = 108 nm, which is 20% 
smaller than [(Z£) _1 + (l*)^ 1 = 131 nm. 



B. PHONON OCCUPATION 

In Fig. [SI we show an example of the computed phonon 
occupation n(k,x) at voltage 1.4 V, £=300 nm. The 
phonon occupation is clearly distributed in a not uniform 
way. Two important peaks with n ~ 20 are present near 
the electrodes. For a given phonon occupation n(k,x), 
we define an effective temperature T e ^(k,x) as the cor- 
responding temperature within the Bose-Einstein statis- 
tics: n — l/{exp[fitu/(kBT e ff)] — 1}. n ~ 20 is extremely 
high and corresponds to T e " ~ 40000 K. Using simple 
arguments, it can be shown that n(fc, x) ^ only for 
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FIG. 5: (Color online) Phonon occupation n as a function 
of the position along the tube, x, and momentum k. Upper 
panel: r-E2 9 LO phonon. Lower panel: K.-A[. The tube is 
the same as in Fig.|21 with voltage 1.4 V. 



fc| < 2Ve/(hvp). This region is small compared to the 
Brillouin zone of graphite 23 . Thus, the hot phonons are 
indeed concentrated around T and K. To evaluate the 
impact of the large phonon occupation on the nanotube 
stability we compute the average kinetic energy (per car- 
bon atom) associated the hot-phonons as 



E = ^ 



v hw v 



2nN at 



n"(k, x) dk dx, 



where N at is the number of atoms in the tube, mF = 1 
and m K = 2 (to take into account both K and K'). 
E at /k B = 50.6 K, at bias 1.4 V and L=300 nm. Thus, 
the actual average kinetic energy per atom, associated to 
the hot phonons, is not particularly high despite the high 
values reached by the hot phonons distribution, n"(k,x). 
Indeed, only a small subset of optical phonons is heated in 
a limited region of the k-space. Concluding, the presence 
of hot phonons having temperatures of several thousands 
K is compatible with the fact that the tubes are not 
melting during electron transport. 

As suggested in Ref. 8, the occupation of the Fi2 g 
phonon at A; = 0, n T (0,x), is accessible experimentally 
using Raman spectroscopy. Indeed, this mode is respon- 
sible for the G peak observed in Raman spectroscopy^ 2 ^. 
The Stokes and anti-Stokes G peak intensities are pro- 
portional to n r (0, x) + 1 and n r (0, x), respectively. From 
the measured Stokes and anti-Stokes ratio one can ex- 



tract n r (Q,x) of a SWNT during electron transport. In 
the upper panel of Fig. [21 we show, as a function of the 
applied bias, the effective temperature T e ^ of the Ra- 
man active E2 9 phonon, obtained from the average over 
the tube length of n r (0,x). Already at 1.0 V, we find a 
T e H = 2500 K (Fig. corresponding to n = 0.7. 

Finally, Refs. l8l9l built two oversimplified models to de- 
scribe the hot-phonon generation. Both models are based 
on the assumption that there is a simple relation between 
the effective scattering length l sc , the zero-temperature 
back-scattering length lb s , and the average phonon occu- 
pation n: 



2n + l 



(6) 



(see Eq.10 of Ref.H and lines 4-12, page 2 of Ref. 0). Al- 
htough a priori reasonable, Eq.[B]is in contradiction with 
the observation that the differential resistance (and hence 
l sc ) is constant with the voltage at high-bias (Fig. [3}, 
while the phonon occupation n is increasing with the 
voltage (upper panel of Fig.[2J). Therefore Eq. IBJcan not 
be used to construct quantitative models for the descrip- 
tion of hot phonons in metallic SWCNs. The high-bias 
linearity of the IV curve, that is observed experi menta lly 
in tubes with L < 500 nm, might have led Refs. Irld lTllOl 
into neglecting hot-phonons effects for tubes lying on a 
substrate. 



IV. CONCLUSIONS 

In conclusion, the explicit and self-consistent dynamics 
of both of phonons and electrons is the necessary ingredi- 
ent for the reliable description of the hot phonons gener- 
ation during electron transport in metallic SWNTs. We 
find a remarkable agreement with the experimental IV 
curves. We have shown that the hot phonon generation 
can not be properly described using the simplified mod- 
els reported in Refs. We demonstrate the presence of 
hot phonons for voltages > 0.2 V also for tubes which are 
lying on a substrate. We have shown that the high-bias 
resistance strongly depends on the thermalization time of 
the optical phonons. As a direct consequence, an impor- 
tant improvement of metallic nanotubes performances is, 
in principle, possible by increasing the optical-phonons 
thermalization time. 
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